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Proton acceleration by using a 620-TW, 18-J laser pulse of peak intensity of 5 x 10^^ W/cm^ 
irradiating a disk target is examined using three-dimensional particle-in-cell simulations. It is shown 
that protons are accelerated efficiently to high energy for a "light" material in the first layer of 
a double- layer target, because a strongly inhomogeneous expansion of the first layer occurs by 
a Coulomb explosion within such a material. Moreover, a large movement of the first layer for 
the accelerated protons is produced by radiation-pressure-dominant acceleration. A time-varying 

CO electric potential produced by this expanding and moving ion cloud accelerates protons effectively. 

^"^ In addition, using the best material for the target, one can generate a proton beam with an energy 

^^ of 200 MeV and an energy spread of 2%. 

;-H PACS numbers: 52.38.Kd, 29.25.Ni, 52.65.Rr 
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INTRODUCTION 



Recently, there has been great progress in compact laser systems, with dramatic improvements in both laser power 

and peak intensity. Ion acceleration by laser pulses has proved to be very useful in applications using compact laser 

O systems. Laser-driven fast ions are expected to be useful in many applications such as hadron therapy,^ ^ fast ignition 

I for thermonuclear fusion,^ ^ laser-driven heavy ion colliders,^ ^ and other applications that use the high-energy ions. 

H Although the achieved proton energy at present is not high enough for some applications such as hadron therapy, 

c/^ which requires 200-MeV protons, other methods can be considered for generating higher energy protons. One simple 

^ ^ way is by using a higher power laser. However, current power capabilities of compact lasers are insufficient; moreover, 

f^ laser power enhancement will result in a cost increase of the accelerator. Therefore, it is important to study conditions 

(y5 for generating higher energy protons with lower laser power and energy by using some special techniques .^^Hi^ 

O In this paper, I show a way to obtain 200-MeV protons by using a laser pulse whose intensity is Iq ~ 10^^ W/cm^, 

*(y^ energy is Eias < 20 J, and power is P ^ 500 TW. I use three-dimensional (3D) particle-in-cell (PIC) simulations to 

K^ investigate how high-energy, high-quality protons can be generated by a several- hundred- terawatt laser. I study the 

F^ proton acceleration during the interaction of the laser pulse with a double-layer target composed of a high-Z atom 

, ^~ ^ layer coated with a hydrogen layer (see Fig. fl]). As suggested in Refs. T and TBHIOl a quasimonoenergetic ion beam 

can be obtained using targets of this type. Our aim is to obtain a high-energy {£ ^ 200 MeV) and high-quality 

^-H {A£/S < 2%) proton beam using a relatively moderate power laser. 

^ In the following sections, I show the dependence of the proton energy on the material of the first layer and that the 

^n high-energy protons can be generated by optimally combining a couple of ion acceleration schemes. 
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II. ION ACCELERATION 

I consider ion acceleration by a charged disk. The charged disk is produced by a laser pulse with sufficiently high 
intensity irradiating a thin foil. Many electrons are driven from the foil by the laser pulse, although the ions of the 
foils almost stay at their initial positions because they are much heavier than the electrons. Therefore, the thin foil 
will have a charge, which induces an electrostatic field. Ions located on the foil surface are accelerated by this electric 
field. The x component of the electric field of a positively charged thin disk is 

where p is the charge density, / is the disk thickness, eo is the vacuum permittivity, and R is the charged disk radius. 
I assume that the x axis is normal to the disk surface placed at the disk center. The solid curve in Fig. [l] shows this 
electric field. The ions, i.e. protons, are accelerated in this electric field, although it rapidly decreases as a function 
of distance from the target surface. The electric field decreases to 10% at x = 2i?, which is the distance equal to the 
diameter of the target and can be considered to be the spot size of a laser pulse. Therefore, generating higher energy 
protons requires producing a higher surface charge density, p/, or increasing R. The former requires a higher intensity 
laser and the latter requires a higher power laser. The rapidly decreasing accelerating field and its narrow width lead 




FIG. 1: Configuration of a double-layer target. The x component of the electric field, Ex(x), is normalized by its maximum, 
pi /2eo, of an electrically charged disk on the x axis (solid curve). Protons are accelerated in this electric field. 

to inefficient proton acceleration by the charged disk. In this paper, I present ways to improve these inefficiences and 
to generate high-energy protons effectively. 

Here, let us define the some terms. In laser ion acceleration, the ions are accelerated in some electric field, E. We 
assume that for an ion of mass of m and charge q^ the force on it from the electric field is qE. The equation of motion 
is qE = -^{mv), where v is the ion velocity. This equation can be written as 



E 



dt 



{mv)^ 



(2) 



where m = m/q. m is the resistance to movement of an ion in a certain electric field, E; therefore we call m "mass" 
in this paper. This expression shows that the smaller rh ions can experience greater acceleration in a certain electric 
field, E. Therefore, small- "mass" ions will be called "light," and big- "mass" ions will be called "heavy." Ions of 
the same "mass" undergo the same movement in a certain electric field. Note that m is equal to the inverse of the 
well-known parameter q/m^ the charge-to-mass ratio. I use m in this paper because it makes it very simple and easy 
to image the movement of charged particles in an electric field. 

Figure [T] shows that the accelerating protons exit the electric field in a short time. This means that the electric field 
produced is not used enough for proton acceleration. Therefore, we should create a situation in which the protons 
experience this electric field longer for efficient acceleration. If the electric potential moves in the direction of the 
moving protons, the protons will experience the electric field longer. In other words, the charged first layer keeps 
pushing the moving protons. I present two ways to create this situation. 

One way this situation can be created by the use of a Coulomb explosion of the first layer. Figure |2] shows that 
the first layer disk undergoes a strongly inhomogeneous expansion owing to the Coulomb explosion. This expansion 
raises the moving electric potential for the accelerating protons. In other words, many ions in the first layer are 
distributed close to the accelerating protons keeping a comparatively high density and move in the proton direction. 
The acceleration rate is higher when the expansion velocity of the first layer ions is higher. This means that the 
strong Coulomb expansion operates effectively for proton acceleration. The Coulomb explosion level is determined by 
the "mass" of the ions composing the first layer. Equation (pi shows that "light" ions have a high expansion velocity. 
That is, "light" ions undergo a stronger a Coulomb explosion and should be generating higher energy protons. 

Another way to induce movement of the first layer is by radiation-pressure-dominant acceleration (RPDA). Figure 
p] shows that the first layer, which expands by a Coulomb explosion (with an ellipsoidal light pattern), is moving with 
velocity V in the laser propagation direction (proton direction) by RPDA. This movement leads the moving electric 
potential. Higher V values generate higher energy protons, since the protons experience the accelerating electric field 
over a longer time by following the electric potential. A portion of the energy and momentum transferred from the 
laser pulse to the electrons is imparted to the ions via a charge separation field. That is, the ions get accelerated by 
this field, and the "light" ions have higher velocity (Eq. ([2])). Thus the "light" ions experience a higher first layer 
velocity and should be generating higher energy protons. 



^ 




Exploding 
first layer 




Accelerating 
protons 



Moving 
electric 



Gentler of 
1 \;the exploding 
'f first layer 



potential 



Initial t 
first layer 



X 



ui 



FIG. 2: The first layer using "ligiit" materials produces a strongly inhomogeneous expansion due to the Coulomb explosion 
(light pattern). The expanding first layer moves at average velocity V in the direction of laser propagation by RPDA. The 
electric potential moves in the x direction as a result of these effects. 

One can obtain higher energy protons by using a "light" material in the first layer, as is corroborated by the 
simulations described below. The simulations were performed with a 3D massively parallel electromagnetic code, 
based on the PIC method.^ 

III. SIMULATION OF FIRST LAYER MATERIALS 

In this section, I study the dependence of the proton energy on the first layer materials of the double-layer target 
by using simulations. 

A. Simulation parameters 



Here, I show the parameters used in the simulations. The spatial coordinates are normalized by the laser wavelength 
A = 0.8 /im and time is measured in terms of the laser period, 27r /uj. 

I use an idealized model, in which a Gaussian linearly polarized laser pulse is incident on a double-layer target 
represented by a collisionless plasma. 

The laser pulse with dimensionless amplitude a = qeEo/rrieUJC = 50, which corresponds to a laser peak intensity 
of 5 X 10^^ W/cm^, is lOA long in the propagation direction, 27 fs in duration, and focused to a spot with size 4A 
(FWHM), which corresponds to a laser peak power of 620 TW and a laser energy of 18 J. Here, the laser peak power 
is calculated by using /^ J^ /(?/, z)dydz and the laser energy is calculated by using J^ J^ J^ /(?/, z, t)dydzdt, 
where / is the laser intensity. The laser pulse is normally incident on the target. The electric field is oriented in the 
y direction. I use this laser pulse in all simulations in this paper. 

Both layers of the double-layer target are shaped as disks. The first layer has a diameter of 8A and a thickness of 
0.5A. The second, hydrogen, layer is narrower and thinner; its diameter is 4A and its thickness is 0.03A. The electron 
density inside the first layer is rig = 3 x 10^^ cm~^ and inside the hydrogen layer it is ng = 9 x 10^^ cm~^. The 
total number of quasiparticles is 8 x 10^. The number of grid cells is equal to 3300 x 1024 x 1024 along the X, y, 
and Z axes, respectively. Correspondingly, the simulation box size is 120 A x 36. 5 A x 36. 5 A. The boundary conditions 
for the particles and for the fields are periodic in the transverse (Y^Z) directions and absorbing at the boundaries of 
the computation box along the X axis, xyz coordinates are used in the text and figures; the origin of the coordinate 
system is located at the center of the rear surface of the initial first layer, and the directions of the x, ?/, and z axes 
are the same as those of the X, y, and Z axes, respectively. That is, the x axis denotes the direction perpendicular 
to the target surface and the y and z axes lie in the plane of the target surface. 

Although the first layer material can be varied, the number of ions is the same in all cases and the ionization state 
of each ion is assumed to be Zi = -\-6 ^^^ since the laser parameters and the target geometry are fixed. 
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FIG. 3: Particle distribution and electric field magnitude (isosurface for value a = 2) for carbon (a) and gold (b). Half of 
the electric field box has been removed to reveal the internal structure. Shown are the initial shape of the target and the laser 
pulse (t = 0), the interaction of the target and laser pulse (t = 25,50 x 27t/uj), and the first layer shape and the accelerated 
protons (color scale) (t = 75, 100 x 27r/cj). 



B. Proton energy as a function of first layer materials 



To examine the dependence of the proton energy, £^ on the first layer material, I performed simulations using 
different materials at normal incidence to the laser pulse. First, I show two simulation results, one in which the first 
layer consists of carbon (a "light" material) and the other is in which it consists of gold (a "heavy" material). The 
results are shown in Figs. |3}|8] 

Figure [3] shows the particle distribution and electric field magnitude of both cases at each time. We see that the 
carbon ions are distributed over a much wider area by its Coulomb explosion than are the gold ions. However, the 
deformation of the laser pulse is similar in both cases at all times. The Coulomb explosion process of the first layer 
is much slower than the laser pulse progress. The first layer is almost undeformed at t = 25 x 27t/uj when the laser 
pulse is just around the target and it has the strong interactions with the target. A big deformation of the first layer 
appears at t > 25 x 27t/uj when the laser pulse passes through or reflects from the target. Therefore, the explosion of 
the first layer is an almost simple Coulomb explosion without other effects. The proton energy obtained in the carbon 
case is much higher than that in the gold case. The average proton energy at t = 100 x 2i\ juj^ £^ave, is 63 MeV for the 
carbon case and 38 MeV for the gold case (a factor of 1.7 higher). 

The reason why higher energy protons are obtained in the carbon case is that the first layer deformation effectively 
contributes to the proton acceleration. Let us examine the deformation of the first layer. Figure [4] shows a cross 
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FIG. 4: Distribution of the ions of the first layer and protons (color scale) in the carbon (a) and gold (b) cases at t = 100 x 27r/cj; 
a two-dimensional projection is shown looking along the z axis. In the carbon case, the distribution area of the carbon ion 
cloud is very wide and it moves in the x direction. 



section of the ion density distribution near the (x, y, z = 0) plane at t = 100 x 27r/co', as seen by looking along the z 
axis. In the carbon case, the strong Coulomb explosion distributes the carbon ions ellipsoidally, and the distribution 
area is much wider than in the gold case. The expansion of the cloud of carbon ions is strongly inhomogeneous and it 
appears to be substantially elongated in the longitudinal direction. The surface of the carbon ion cloud is close to the 
acceleration protons, which means that the protons keep getting pushed by a comparatively strong force. Moreover, 
the center point of the ion cloud is moving in the x direction, at a coordinate of 3.1 A, so that the electrostatic potential 
is moving in the direction of the protons. In contrast, in the gold case, the ion distribution is very compact and the 
distance between the ion cloud surface and the protons is much greater than in the carbon case. In addition, it almost 
does not move, with the coordinate of the center point of gold ions being O.OA. Moreover, the proton positions in the 
carbon case have travelled further along the x axis than in the gold case, since the protons in the carbon case have 
higher energy (higher velocity) than in the gold case. 

Here, let us consider the variation of the proton energy in time for both cases. Figure Isja) shows the average 
proton energy versus time. The proton energy in the carbon case is always higher than in the gold case and difference 
in energy grows in time. The protons are accelerated relatively quickly until a time of about t = 50 x 2i\ juj in 
both cases. In the gold case, the acceleration almost saturates at t > 50 x 2i\ juj^ although in the carbon case it is 
still increasing compared with the gold case. This is because the electric field seen by protons in the carbon case is 
greater and continues longer than in the gold case, owing to the expansion and movement of first layer (see Fig. [4|. 
Correspondingly, their energy gain is also greater. In the carbon case the proton energy is 1.7 times higher than in 
the gold case and the energy spread is almost the same in the both cases, as seen in Fig. [5|b). The energy spread, 
Af /fave, is 6% and 8% for the cases of carbon and gold, respectively. Higher energy protons can be obtained in the 
carbon first layer, as a result of the stronger electric field and movement of the first layer. I consider the differences 
in the electric field and the first layer movement in the two cases. 

First, I show the differences in the electric field between the two cases. Figure [6] shows the x component of the 
electric field along the x coordinate (solid lines), Ex(x\ and the proton position (circular dots) at each time. In the 
gold case (see Fig. [6]^b)), the electric field simply decreases with the distance from the surface of the initial target at 
all times. Therefore, the electric field that protons experience (see the circular dots) simply decreases too. In contrast. 
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FIG. 5: (a) Average proton energy S versus time, as obtained in the simulation shown in Figs. Is] The proton energy 
of the carbon case is higher than in the gold case at all times, (b) Proton energy spectrum obtained in the simulation at 
t — 100 X 2ti /uj^ for the cases of gold and carbon, normalized by the maximum in the former case. 
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FIG. 6: The x component of the electric field, Ex{x) (solid lines) and the proton position (circles) at each time. In the gold 
case, in the inset, the electric field simply decreases with distance from the target surface. In the carbon case, it does not 
decrease much and has a peak that moves to the right side. 



in the carbon case (Fig. [oFa)), the electric field peaks at points near the protons for each time. The electric field 
moves in the direction of the accelerating protons keeping the same shape (see Fig. [6Fa): t = 7b and t — 100 x 27r/cj). 
The electric fields at the proton positions slowly decrease, and those are higher than in the gold case at all times. At 
the proton position at t = 50 x 2i\ juj^ the value of E^ is 2.4 and 1.4 TV/m for the cases of carbon and gold ions, 
respectively. 

Next, I consider the movement of the first layer. Figure [7| shows the laser pulse and the target at an early time, 
t < 20 X 27r/cc;, for the carbon case, as seen by looking along the z axis. The laser pulse moves from the left-hand 
size to the right-hand side, and the color shows the electric field magnitude. We can see that many electrons (gas 
pattern) are pushed out from the target by the ponderomotive force from the laser pulse. Those pushed-out electrons 
are distributed to the rear area (propagation direction) of the target and they move forward in the laser propagation 
direction, although the carbon ions stay at their initial position. This charge separation produces the strong electric 
field. Then, the carbon ions are moved in the laser propagation direction by this electric field (which is RPDA). 
Figure [8] shows the first layer velocity, V^^ normalized by the speed of light, and the first layer position, normalized 
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FIG. 7: Spatial distribution of particles and the electric field magnitude at early simulation times for the carbon case. A 
two-dimensional projection is shown looking along the z axis. The solid bar is a projection of the first layer (t = 0). The gas 
pattern denotes electrons that are pushed out from the target by the laser pulse (t = 10, 20 x 2t\ juS) and distributed to the rear 
area (propagation direction) of the target. This produces the movement of the first layer ions. 
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FIG. 8: (a) Velocity of the first layer of the target in the x direction normalized by the speed of light, \4/c, as a function of 
time. The velocity of the carbon target is much higher than that for the gold case, (b) Movements of the first layer in the x 
direction normalized by the wavelength, x/A, as a function of time. 



by the wavelength, for the x direction as a function of time. These are averaged values for all ions of the first layer. 
The first layer velocity rises rapidly at the initial time, t ~ 20 x 2i\ juj^ when the laser pulse is still around the target 
(see Fig. ItI) and the velocity is constant at time t > 25 x 27r/co', after the laser pulse passes through or reflects off the 
target. The increase in the fist layer velocity stops at t ^ 20 x 27r/cj, since by this time there is no clearly one-sided 
distribution of the electrons like that at time t ~ 10 x 2i\ juj (see Fig. ItI). The first layer velocity in the carbon case is 
14 times that in the gold case at t > 25 x 2i\ juj. This means strong RPDA occurs in the carbon case. The movement 
in the carbon case is much greater than in the gold case too, and the difference in distance between the carbon case 
and the gold case grows with time. This indicates that the moving electric potential operates efficiently in the carbon 
case. Incidentally, the velocity for the y direction, Vy^ of the first layer is relatively very small. It is about 1/1000 of 
Vrj. in both the carbon and gold cases. Movement for the y direction is very small too, because they are normal to the 
incidence of the laser pulse. 

In above considerations, I showed that the higher energy protons can be obtained for "lighter" material in the first 
layer by comparing carbon and gold. Here, I investigate the effect using additional materials. Figure [9] shows the 
average proton energy for different materials comprising the first layer. The proton energy is normalized by the energy 
in the case of gold ions. We see that the average proton energy almost linearly depends on the ratio of Z\l A^ where 
Z\ is the charge state and A is the atomic mass number of the ion. Higher energy protons can be obtained by using 
a larger ratio of Z^/A^ i.e., "light," material for the first layer in the double-layer target. 

In experiments on laser-driven ion acceleration, a CH polymer target exhibited higher energy protons^^ than a 
metallic target.— 
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FIG. 9: Average proton energy, normalized by energy in the gold case, shown with respect to the ratio between the charge 
state, Zi, and the atomic number, A, of ions making up the first layer. 



IV. COMPONENTS OF THE PROTON ENERGY 



In the previous section, I showed that higher energy protons can be obtained by using "light" material in the first 
layer, because a strong Coulomb explosion and a greater first layer movement toward the accelerating protons occur 
in such materials. These effects augment the acceleration by the electric field of the charged first layer disk. In this 
section, I show the amount of each effect of the acceleration on the total proton energy in the carbon case. 

To examine the effect of the first layer velocity, V ^ I performed simulations with a laser pulse by reversing the 
irradiation direction of the laser (reverse irradiating) in the previous carbon case. Thus the hydrogen layer was put 
on the front side of the target (see Fig. 10). The results are shown in Figs. 10 and 11 The protons are accelerating 



in the +x direction even when the laser pulse irradiates the target in the reversed way. The average proton energy 
in this case, £~ ^ is 32 MeV. Since the proton layer is very thin and small, it has less effect on the first layer velocity. 
Therefore, the first layer velocity, V~ ^ has the same absolute value as in positive irradiation, the case of the previous 
section, but the sign is opposite, V~{t) — —V{t)^ where V{t) is the first layer velocity in time of the positive case. 
Therefore, the first layer movement has a negative effect on the proton energy in the reverse irradiating case. Using 
this result, we can estimate the amount of proton energy attributable to the first layer velocity, £y. The proton 
energy in the positive irradiating case, f +, and in the reverse irradiation case, £~ ^ are written as 



f + = £:« + Sv 



(3) 



£- 



£a — £\ 



(4) 



where 8a is the proton energy without the work done by the first layer velocity (i.e., attributable to the acceleration 
by the electric field of the charged disk of the first layer and the Coulomb explosion of the first layer ions). Hence we 
obtain the work done by the first layer velocity, Sy = {£^ ~^~)/'^-> which is 16 MeV. Therefore, the ratio of the effect 
of the first layer velocity, RPDA, is 25% of the total proton energy in the carbon case. In our simulations, the laser 
intensity and energy of /q = 5 x lO^^W/cm^ and Eias = 18 J, respectively, are not enough for the RPDA regime in 
full scale, but this shows that the RPDA regime has a strong effect even at this laser power level. The proton energy 
without the work done by the first layer velocity, 8ai is 48 MeV. 

Next, I examine the work done by the Coulomb explosion of the first layer. To do so, I performed a simulation 
with a laser pulse reverse irradiating in the gold case. In this case, the average proton energy, £^ ~, is 32 MeV. In 
the positive irradiation case, the gold case of the previous section, the average proton energy, £ +, was 38 MeV. The 
difference between these values is very small, because the first layer velocity, F, is very small (see Figs. |4]and[8|. 
This means that the work done by the first layer velocity is negligibly small in the gold case. The Coulomb explosion 
effect is negligible, as shown in previous section, too. Therefore, we estimate that the work without the first layer 
velocity and Coulomb explosion in the carbon case, £q ^ {8 ^ ^ £ ~)/2, is 35 MeV. This is almost the energy from 
acceleration by only the electric field of the charged disk of the first layer. The ratio of this effect in the total energy 
of the carbon case is 55%. Then we obtain the work done by the Coulomb explosion as ^^c = £^ci — ^o = 13 MeV. The 
ratio of the Coulomb explosion effect, £c /£^ ^ is 20% in the carbon case (see Fig. 12). 
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FIG. 10: (a) 3D view and (b) the projection onto the {x,y) plane for a laser pulse irradiating onto the other side, the proton 
layer side, of a double-layer target. Shown are the particle distribution and electric field magnitude (isosurface for value a = 2); 
half of the electric field box has been removed. The protons are accelerated in the -\-x direction, and the C ions move in the 
—X direction by RPDA. 
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FIG. 11: Proton energy spectrum obtained in the simulation at t = 100 x 2i^ juj^ for the reverse irradiation case, C 
positive irradiation (previous carbon) case, C^, normalized by the maximum in the former case. 
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FIG. 12: Proton energy of each acceleration mechanism in the carbon case. The proton energy produced by the electric field 
of the nonmoving first layer is almost half the total proton energy, and RPDA and the Coulomb explosion each amount to 
almost a quarter of the total proton energy. 



Another effect to consider is that of the protons being dragged by the electrons that are pushed out of the target 
by the laser pulse. I estimate this effect by using the gold case results. The electrons are mainly pushed out in the 
laser propagation direction. Therefore, in the reverse irradiation case, the dragging effect by electrons acts to reduce 
the proton energy. The proton energy in the positive irradiating case, E +, and in the reverse irradiation case, E ~, 
are written as £^ + = E^ -\- Ey -\- E^ and E ~ = E^ — Ey — E^^ where Ey is the work done by the first layer velocity, 
E^ is the work done by the dragging by the electrons, and E^ is other work. The effect of the first layer velocity 
plus the dragging by electrons, Ey -\- E^ d = {E ~^ — E ~) /2^ is 3 MeV. Since Ey > 0, the effect of the dragging by the 
electrons is E^ < 3 MeV, and the ratio is less than 5% of the total energy of the carbon case. Therefore, the effect 
whereby the protons are dragged by the electrons is very small in our simulations. This is because the electrons move 
very fast compared with the ions and protons, the distance between electrons and protons quickly becomes very large 
compared with that between the first layer ions and protons, and, moreover, the electrons are distributed over a very 
wide area. The electric force decreases with distance by second order. Therefore, the force the electrons exert on the 
protons is very small compared with the force exerted by the ions. 

Next, I consider the theory for the work done by the first layer velocity. I assume that the proton velocity v > V 
and v'^ <C c^. The theoretical formula^ is written as Ey = mV'^{l + ^/2EoJmV^^^) ^ where m is the proton mass 
and Eq is the proton energy in the case of a nonmoving first layer. Denoting the proton velocity in the case of a 
nonmoving first layer by vq and assuming Vq <C c^ we obtain 



£v = 2V{V^Vl^V^), 



(5) 



where Ey is the normalized proton energ y of the work done by the first layer velocity, Ev/Eq^ and V is the normalized 
velocity of the first layer, V/vq. In Fig. 13, I present this theoretical dependence of Ey on V", by using formula ([5). 
The proton energy by RPDA, Ey, grows rapidly with increasing first layer velocity V. The simulation results (circle 
and square dot) are plotted on the figure. The theoretical calculations and the simulation results agree well. Since 
the first layer velocity in the carbon case is not so high, V/vq = 0.13, Fig. 13 shows that if we produce a high velocity 
in the first layer, protons with a few times higher energy can be obtained. In addition, assuming F^ <C 1 we obtain 
the simpler formula Ey = 2^. 
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FIG. 13: Proton energy produced by RPDA, 8v ^ as a function of first layer velocity, V . The proton energy, Sv ^ is normalized 
by the proton energy of the nonmoving first layer case, 8v /^q^ and the first layer velocity, V, is normalized by the proton 
velocity of the nonmoving first layer case, V/vq. The theoretical result (solid line) is given by equation ([5]). The simulation 
results for the gold and carbon cases are plotted with a square and a circle, respectively. 



GENERATING A 200-MEV PROTON BEAM 



I show the way to obtain a 200-MeV proton beam, using the same laser pulse as in the previous section with the 
same normal incidence. In the previous section, I showed that higher energy protons can be obtained by using "light" 
material in the first layer. The "lightest" material is hydrogen. Therefore, we could get higher energy protons by using 
hydrogen for the first layer. I evaluate this contribution by simulation. In this case, even if we put the second layer of 
thin hydrogen on the hydrogen first layer, the second layer has no meaning. Therefore, I use a simple hydrogen disk, 
without a second layer. 

The hydrogen disk target has the same shape and size as that of the first layer of the double-layer target used 
previously. The electron density inside the target is ng = 9 x 10^^ cm~^. Because the hydrogen cloud generated 
by the target is distributed over a wider area than in the previous case, I define a wider simulation box for Y and 
Z directions. The number of grid cells is equal to 5000 x 3000 x 3000 along the X, F, and Z axes, respectively. 
Correspondingly, the simulation box size is 113 A x 67. 5 A x 67. 5 A. The total number of quasiparticles is 4 x 10^. The 
other simulation parameters are the same as those used in previous sections. 

Figure p^a) shows the particle distribution and the electric field magnitude in time. At t = 25 x 27r/cj, the laser 
pulse is just around the target and it has the strong interactions with a target. The target maintains its initial disk 
shape at this time. After t = 25 x 27r/co', the laser pulse passes through or reflects off of the target, and the proton cloud 
produced by Coulomb explosion is growing in time. Figure fl4[b) shows a cross section of the ion cloud at each time. 
Hydrogen ions (protons) are classified by color in terms of energy. We see that the exploded hydrogen disk (hydrogen 
ions) is distributed over a very wide area. The expansion of the cloud of hydrogen ions is inhomogeneous and the 
cloud is elongated in the longitudinal direction. Figure 15 shows the proton energy spectrum at t = 100 x 2i\ juj. The 
maximum energy is £^max = 200 MeV and the average energy is £^ave = 25 MeV. The vertical axis is given in units of 
number of protons per 1 MeV. We can estimate the number of obtained protons by using the required energy and 
energy width from Fig. [15] We see that we can obtain five times higher proton energy than in the gold case, even by 
using the same laser pulse, by using an optimum target material. 

Figure [16] shows the target average velocity, normalized by the speed of light, and the target average position, 
normalized by the wavelength, in the x direction as a function of time. The target velocity rises rapidly at the initial 
time, t ~ 20 X 27r/co', when the laser pulse is still around the target and the velocity is constant at time t > 25 x 27r/a;, 
after the laser pulse passes through or reflects off of the target. This is similar to the carbon case (Fig. [8|, although 
the target velocity of this case is 2.5 times that in the carbon case at time t > 25 x 2tt /uj. That means strong RPDA 
appears in this case, because the target velocity is attributable to RPDA. The movement of the target is much greater 
than in the carbon case too. In the hydrogen disk target, the protons are accelerated more efficiently than in the 
carbon case, because hydrogen is much "lighter" than carbon. 

Because I used the same laser pulse in all simulations, the numbers of electrons pushed out from the target are 
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FIG. 14: (a) Particle distribution and electric field magnitude (isosurface for value a = 2), showing the initial shape of the 
target and the laser pulse (t = 0) and the interaction of the target and laser pulse (t = 25, 50 x 27r/cj). Half of the electric 
field box has been removed to reveal the internal structure. For protons, the color corresponds to energy, (b) Two-dimensional 
projection of the particle distribution, shown looking along the z axis. Half of the proton cloud has been removed to reveal the 
internal structure. 
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FIG. 15: Proton energy spectrum obtained in the simulation at t = 100 x 27t/uj. 
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FIG. 16: (a) Velocity of the hydrogen target in the x direction, normahzed by the speed of Ught, Vx/c, as a function of time, 
(b) Movement of the hydrogen target in the x direction, normahzed by the wavelength, x/A, as a function of time. 
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FIG. 17: Distributions of protons at t = 100 x 27r/(jU and the selected area of the proton bunch. 



estimated to be almost the same in carbon and hydrogen disk cases at the initial simulation time. Therefore, the 
target surface charge is almost the same in the two cases, and the proton energy by the charged disk electric field can 



be estimated using the same considerations as in Section IV, yielding 35 MeV. Therefore, we can estimate that the 
proton energy by RPDA and Coulomb explosion is ^165 MeV. We can say that the RPDA and Coulomb explosion 
effect is much stronger in the hydrogen disk case compared with the carbon case. 



We see that protons are distributed in different areas based on each energy level (see Fig. 14). Moreover, high-energy 



protons are distributed on the -\-x side edge of the proton cloud and are moving in the -\-x direction. Therefore, we 
could select only high-energy protons by using a pinhole with a shutter. The shutter must have high enough accuracy 
and the ability to shield unwanted particles and radiation. The shutter speed must be very fast near the target. The 
accuracy becomes coarser if the shutter position moves to a position farther from the target, because the distance 
between the high-energy protons and the low-energy protons grows in time. It may be difficult to construct a shutter 
that satisfies both the timing accuracy and shielding ability, however, we could get a similar result by using a magnet 
and a slit. The path of a proton is changed for each energy level by a magnetic field, and high-energy protons can be 
taken out by passing through a slit. I show the results using the previous method. The cutoff position is shown in Fig. 
17 Figure [is] shows energy spectrum of cutoff protons at t = 100 x 27r/uj. We obtain a proton beam with a maximum 
energy of fmax = 200 MeV and an average energy of £^ave = 193 MeV with an energy spread of Af /£^ave = 2.3% and 
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FIG. 18: Proton energy spectrum obtained by cutting off the proton bunch at t = 100 x 27r/cj, normaUzed by the raaxiraura. 

a particle number of 2.6 x 10^. This proton beam has high enough energy and quahty for some apphcations (e.g., in 
medical applications). 

VI. CONCLUSIONS 

Proton acceleration driven by a laser pulse irradiating a disk target is investigated with the help of 3D PIC 
simulations. I have found higher energy protons are obtained by using "light" materials for the target. As seen in 
simulations, for these materials a strongly inhomogeneous expansion of the disk target occurs owing to the Coulomb 
explosion, which plays an important role, and RPDA has a strong effect. The time-varying electric potential of the 
inhomogeneous expanding ion cloud and the movement of the ion cloud for the protons efficiently accelerate protons. 
The proton beam energy can be substantially increased by using a "light" material for the target. In our simulations, 
the laser intensity and energy, Iq = 5 x 10^^ W/cm^ and Sias = 18 J, are not enough to reach the RPDA regime in 
full scale, but the RPDA regime has a big effect even at this laser power level. Although I show simulation results 
by using a simple hydrogen disk, it may be difficult to fabricate such a hydrogen disk; a CH^ foil target with a high 
n value should be a good substitute. The laser parameters used in this paper — intensity, power, energy, and spot 
size — are ones already existing in current laser systems. Therefore, we should be able to generate 200-MeV protons 
now. 
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